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The Zaremba boundary-value problem is a boundary value prob- 
lem for Laplace-type second-order partial differential operators acting 
on smooth sections of a vector bundle over a smooth compact Rieman- 
nian manifold with smooth boundary but with non-smooth (singular) 
boundary conditions, which include Dirichlet conditions on one part 
of the boundary and Neumann ones on another part of the boundary. 
We study the heat kernel asymptotics of Zaremba boundary value 
problem. The construction of the global parametrix of the heat equa- 
tion is described in detail and the leading parametrix is computed 
explicitly. Some of the first non-trivial coefficients of the heat kernel 
asymptotic expansion are computed explicitly. 
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1 Introduction 

The heat kernel of elliptic partial differential operators acting on sections of 
vector bundles over compact manifolds proved to be of great importance in 
mathematical physics. In particular, the main objects of interest in quantum 
field theory and statistical physics, such as the effective action, the parti- 
tion function, Green functions, and correlation functions, are described by 
the functional determinants and the resolvent of differential operators, which 
can be expressed in terms of the heat kernel. The most important operators 
appearing in physics and geometry are the second order partial differential 
operators of Laplace type; such operators are characterized by a scalar lead- 
ing symbol (even if acting on sections of vector bundles). Within the smooth 
category this problem has been studied extensively during last years (see, for 
example, |30], [10]; for reviews see |9|, (|, |3| and references therein). In the case 



of smooth compact manifolds without boundary the problem of calculation 
of heat kernel asymptotics reduces to a purely computational (algebraic) one 
for which various powerful algorithms have been developed [|l|, ; this prob- 
lem is now well understood. In the case of smooth compact manifolds with 
a smooth boundary and smooth boundary conditions the complexity of the 
problem depends significantly on the type of the boundary conditions. The 
classical smooth boundary problems (Dirichlet, Neumann, or a mixed com- 
bination of those on vector bundles) are the most extensively studied ones 
(see fl2| , |13[ |34] , |2| and the references therein). A more general scheme, so 
called oblique (or Grubb-Gilkey-Smith) boundary value problem |3l], ^D, |28 



which includes tangential (oblique) derivatives along the boundary, has been 
studied in H g H H, II- In this 

case the problem is not automatically 
elliptic; there is a certain strong ellipticity condition on the leading symbol 
of the boundary operator. This problem is much more difficult to handle, 
the main reason being that the heat kernel asymptotics are no longer poly- 
nomial in the jets of the symbols of the differential operator and the bound- 
ary operator. Another class of boundary value problems are characterized 
by essentially non-local boundary conditions, for example, the spectral or 
Atiyah-Patodi-Singer boundary conditions [^0[ [11], |38|| . 

All the above described boundary value problems were smooth. A more 
general (and much more complicated) setting, so called singular boundary 
value problem, arises when either the symbol of the differential operator or the 
symbol of the boundary operator (or the boundary itself) are not smooth. In 
this paper we study a singular boundary value problem for a second order par- 
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tial differential operator of Laplace type when the operator itself has smooth 
coefficients but the boundary operator is not smooth. For the case when the 
manifold as well as the boundary are smooth, but the boundary operator 
jumps from Dirichlet to Neumann on the boundary, is known in the litera- 
ture as Zaremba problem. Such problems often arise in applied mathematics 
and engineering and there are some exact results available for special cases 
(two or three dimensions, specific geometry, etc.) |f42| , |25|j . Zaremba problem 
belongs to a much wider class of singular boundary value problems, i.e. man- 
ifolds with singularities (corners, edges, cones etc.). There is a large body 
of literature on this subject where the problem is studied from an abstract 



function-analytical point of view [26, 17, O, ftM, 14, 15, 41, B7L 35, B3l 27 



However, the study of heat kernel asymptotics of Zaremba type problems is 
quite new, and there are only some preliminary results in this area [f| [21], [20 1. 
Moreover, compared to the smooth category the needed machinery is still un- 
derdeveloped. We would like to stress that we are interested not only in the 
asymptotics of the trace of the heat kernel, i.e. the integrated heat kernel 
diagonal, but also in the local asymptotic expansion of the off-diagonal heat 
kernel. 

This paper is organized as follows. In the sect. 2 a formal description 
of the Zaremba problem is given. In the sect. 3 the general form of the 
heat kernel asymptotic expansion is described. In sect. 4 the construction of 
the global parametrix to the heat equation is described. In sect. 5 the first 
non-trivial heat kernel coefficients are computed explicitly. In conclusion we 
summarize (and discuss) the results and outline some future problems. 



2 General Setup 

2.1 Laplace Type Operators 

Let (M, g) be a smooth compact Riemannian manifold of dimension m with 
a boundary dM, equipped with a positive definite Riemannian metric g. Let 
V be a vector bundle over M, V* be its dual, and End (V) = V <g> V* be 
the corresponding bundle of endomorphisms. Given any vector bundle V, we 
denote by C°°(y) its space of smooth sections. We assume that the vector 
bundle V is equipped with a Hermitian metric. This naturally identifies the 
dual vector bundle V* with V, and defines a natural L 2 inner product and 
the L 2 -trace using the invariant Riemannian measure d vol g on the manifold 
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M. The completion of C°°(V) in this norm defines the Hilbert space L 2 (V) 
of square integrable sections. 

We denote by TM and T*M the tangent and cotangent bundles of M. 
Let a connection, V y : C°°(V) -> C°°(T*M ® 7), on the vector bundle 7 
be given, which we assume to be compatible with the Hermitian metric on 
the vector bundle V. The connection is given its unique natural extension 
to bundles in the tensor algebra over V and V*. In fact, using the Levi- 
Civita connection V LC of the metric g together with V v , we naturally obtain 
connections on all bundles in the tensor algebra over V, V*, TM and T*M; 
the resulting connection will usually be denoted just by V. It is usually 
clear which bundle's connection is being referred to, from the nature of the 
section being acted upon. We also adopt the Einstein convention and sum 
over repeated indices. With our notation, Greek indices, . . ., label the 
local coordinates on M and range from 1 through m, lower case Latin indices 
from the middle of the alphabet, i,j, k,l, . . ., label the local coordinates on 
DM (codimension one manifold) and range from 2 through m, and lower case 
Latin indices from the beginning of the alphabet, a,b,c,d, . . ., label the local 
coordinates on a codimension two manifold S C dM that will be described 
later and range over 3, . . . , m. Further, we will denote by g the induced 
metric on the submanifolds (of the codimension one or two) and by V the 
Levi-Civita connection of the induced metric. We should stress from the 
beginning that we slightly abuse the notation by using the same symbols for 
all submanifolds (of codimension one and two). This should not cause any 
misunderstanding since it is always clear from the context what is meant. 

Let V* be the formal adjoint of the covariant derivative defined using 
the Riemannian metric and the Hermitian structure on V and let Q G 
C°°(End(V)) be a smooth Hermitian section of the endomorphism bun- 
dle End (V). The Laplace type operator F : C°°(V) — > C°°(V) is a partial 
differential operator of the form 

F = V*V + Q = -g' w V li V v + Q. (1) 

Alternatively, the Laplace type operators are second-order partial differ- 
ential operators with positive definite scalar leading symbol of the form 
ctl(F;x, £) = I|£| 2 = lg liv {x)^n^ v . Hereafter I denotes the identity endo- 
morphism of the vector bundle V . We will often omit it whenever it does not 
cause any misunderstanding. Any second order operator with a scalar lead- 
ing symbol can be put in the form (fl|) by choosing the Riemannian metric g, 
the connection V v on the vector bundle V and the endomorphism Q. 
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2.2 Boundary Conditions 

In the case of manifolds with boundary, one has to impose some boundary 
conditions in order to make a (formally self-adjoint) differential operator self- 
adjoint (at least symmetric) and elliptic. Let N be the inward-pointing unit 
normal vector field to the boundary and let W = V\ dM be the restriction of 
the vector bundle V to the boundary dM. We define the boundary data map 
i> : C°°(V) ^L 2 {W®W) by 

The boundary conditions then read 

£#/?) = 0, (3) 

where B : L 2 (W © W) — > L 2 (W © W) is the boundary operator, which will 
be specified later. If the operator B is a tangential differential operator (pos- 
sibly of order zero), then the boundary conditions are local. Otherwise, for 
example, when B is a pseudo-differential operator, the boundary conditions 
are non-local. 

To define the boundary operator one needs a self- adjoint orthogonal pro- 
jector II that splits the space L 2 (W) in two orthogonal subspaces 

L 2 {W) = L\{W) @ L\{W) , (4) 

where 

L 2 (W) =UL 2 {W), and L\ = ( Id - U)L 2 {W) , (5) 

and a self-adjoint operator A : L 2 (W) — > L 2 (W), such that AL 2 ^(W) = 
{0}, i.e. IIA = ALT = 0. Hereafter Id denotes the identity operator. The 
boundary operator is then defined by 

which is equivalent to the following boundary conditions 

:0, (7) 

id - n) (Vjvp )+a(V ) = o, (8) 

V dM/ V dM/ 



n(^ 



dM 
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It is easy to see that the boundary operator B and the operator 

* = M-B=( u _; n £)■ w 

are complimentary projectors on L 2 (W © W 7 ), i.e. 

B 2 = B K 2 = K, BK = KB = . (10) 

Hence, a section that satisfies the boundary conditions can be parametrized 
by x(<P) = u((p) © v((p) G L 2 (iy © iy), u(<p) e i;(<p) G Ljj, so that 

It is not difficult to see that the boundary operator B (|6]) incorporates all 
standard types of boundary conditions. Indeed, by choosing II = I and A = 
one gets the Dirichlet boundary conditions, by choosing II = 0, A = I one 
gets the Neumann boundary conditions. More generally, the choice IT, A G 
C°°( End {W)), so that All = IIA = 0, corresponds to the mixed boundary 
conditions. 

Remark 1 The boundary dM could be, in general, a disconnected manifold 
consisting of a finite number of disjoint connected parts, dM = U" =1 £j ; with 
each Ej being compact connected manifold without boundary, <9£j = and 
£j PI Sj = if i 7^ j. Thus one can impose different boundary conditions 
on different connected parts of the boundary S». This means that the full 
boundary operator decomposes B = B\ © • • • © B n , with Bi being different 
boundary operators acting on different bundles. 

We always assume the manifold M itself and the coefficients of the opera- 
tor F to be smooth in the interior of M. If, in addition, the boundary dM is 
smooth, and the boundary operator B is a differential operator with smooth 
coefficients, then (F, B) is called smooth local boundary value problem. 

In this paper we are interested in a different class of boundary conditions. 
Namely, we do not assume the boundary operator to be smooth. Instead, 
we will study the case when it has discontinuous coefficients. Such prob- 
lems are often called mixed boundary conditions; to avoid misunderstanding 
we will not use this terminology. We impose different boundary conditions 
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on connected parts of the boundary, which makes the boundary value prob- 
lem discontinuous. Roughly speaking, one has a decomposition of a smooth 
boundary in some parts where different types of the boundary conditions are 
imposed, i.e. Dirichlet or Neumann. The boundary operator is then dis- 
continuous at the intersection of these parts. The boundary value problems 
of this type are called Zaremba problem in the literature 15 1 (see also 



In this paper we consider the simplest case when there are just two com- 
ponents. We assume that the boundary of the manifold dM is decomposed 
as the disjoint union 

0M = E 1 UE 2 UE o , (12) 

where Si and E 2 are smooth compact submanifolds of dimension (m — 1) 
(codimension 1 submanifolds), with the same boundary £ = <9£i = d"E 2 , 
that is a smooth compact submanifold of dimension (m — 2) (codimension 2 
submanifold) without boundary , i.e. <9£o = 0. Let us stress here that when 
viewed as sets both Ei and £ 2 are considered to be disjoint open sets, i.e. 

Sxns 2 = 0. 

Let Xi '■ 9M — * R, (i = 0,1,2), be the characteristic functions of the 
sets £j, Xi{&) = 1 if ^ G Sj and Xi(&) = if a; G" £j. Obviously, Xi(%) + 
Xa(z) + Xo(x) = 1 for any x G dM. Let 7T; : L 2 (W) -> L 2 (W), (i = 0, 1,2), 
be the trivial projections of sections, ip, of a vector bundle W to Ej defined 
by (7IY0)(£) = Xi(^)ij(x), i.e. (iriip)(x) = ip(x) if x G Ej and (7r^)(x) = 
if x G" Ej. In other words 7i"i maps smooth sections of the bundle to 
their restriction to Ex, extending them by zero on E 2 , and similarly for 7r 2 . 
Obviously, tti + 7r 2 + tt = Id, tt? = 7Tj, (i = 0, 1, 2), and 7Tj7Tj = for i ^ j. 
Let A G C°°( End (W^)) be a smooth Hermitian endomorphism of the vector 
bundle W. Then the boundary operator of our problem can be written in 
the form 

fl=f 7 ° V (13) 

The projectors tti and 7r 2 as well as the boundary operator B are clearly non- 
smooth {discontinuous) on Eq. In other words, we have Dirichlet boundary 
conditions on E x and Neumann (Robin) boundary conditions on E 2 : 



0, (14) 



(Vtv + A)^ =0 (15) 



7 



Note that the boundary condition are set only on open subsets Si and E 2 ; 
the boundary conditions do not say anything about the boundary data on 
S . We will see later that to specify the solution uniquely we also need a 
further condition which specifies the type of the singularity on Eo- 

We will call the boundary value problem (F, B) for a Laplace operator F 
([]]) with the boundary operator B of the form (|T3| ) Zaremba boundary value 
problem. 

2.3 Symmetry 

Let us define the antisymmetric bilinear form 

I(<Pl, <P2) = (F<Pl, ^2)l?{V) ~ {<Pl, Fip2) L 2(y) , (16) 

for any two smooth sections <fi,if2 G C°°(V) of the vector bundle V. By 
integrating by parts on M one can easily see that this bilinear form depends 
only on the boundary data 

I(<Pl, <P2) = J^M)l^(w®w) , (17) 

where 

(18) 

Therefore, it vanishes on sections of the bundle V with compact support 
disjoint from the boundary dM when the boundary data vanish VK^i) = 
"0(^2) = 0. This is a simple consequence of the fact that the operator F 
is formally self-adjoint. A formally self-adjoint operator is essentially self- 
adjoint if its closure is self-adjoint. This means that the operator is such that: 
i) it is symmetric on smooth sections satisfying the boundary conditions, 
and ii) there exists a unique self-adjoint extension of it. To prove the latter 
property one has to study the deficiency indices; however, this will not be 
the subject of primary interest in the present paper. We check only the first 
property, i.e. that the operator F is symmetric. 

By integrating by parts on dM, it is not difficult to check that the form 
I((fx, <p2) does vanish for any <p\, <p 2 G C°°(V) satisfying the boundary condi- 
tions with the boundary operator (|6|) provided the operator A is symmetric. 
Therefore, we immediately obtain that the Zaremba boundary value problem 
is symmetric. 
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2.4 Ellipticity 

Let F Q be a Laplace type operator with constant coefficients obtained from 
the operator F by freezing the coefficients at a fixed point xo in the leading 
derivative part, i.e. for the Laplace type operator F = —g^ u (xo)d^d u . If 
the point Xq is in the interior of the manifold M, we assume the operator 
F to act on the sections of the vector bundle V over R m . If the point x 
is on the boundary, we assume the operator F to act on sections of V over 
]R m_1 x R + . If the point x is on the boundary, we also define the operator 
i?o by omitting the non-diagonal part of B and freezing the coefficients of 
the diagonal part of the boundary operator B, i.e. B = Xi(^o) © X2 (£())• 
Note that the coefficients of the operator B are constant on Ei and E 2 but 
are discontinuous on £ . 

The boundary value problem (F, B) is said to be elliptic with respect to 
C \ IR+ if for any complex number A G C \ M + not lying on the positive real 
axis the following two conditions are satisfied: i) for any interior point x the 
equation 

(F - A)y? = (19) 

has a unique non-trivial solution vanishing at infinity, and ii) for any bound- 
ary point x the above equation has a unique non-trivial solution vanishing 
at infinity and satisfying the boundary conditions 

B o i/>((p) = 0. (20) 

Here non-trivial solution means that it is not identically zero, tp ^ 0, for 
A 7^ and it is not constant, ip ^ const , in case A = 0. 

The question of ellipticity of Zaremba boundary value problem is a subtle 
one. We will show below that for the freezed problem to have a unique 
solution one has to impose an additional condition along the codimension two 
submanifold E , which specifies the (singular) behavior of he solution near 
So- Without this condition (which is often imposed implicitly by choosing 
the most regular solution) the freezed problem has infinitely many acceptable 
(square integrable) solutions, so that Zaremba problem fails to be elliptic. 

3 Heat Kernel 

For t > the heat semi-group operator U(t) = exp(-tF) : L 2 (V, M) — > 
L 2 (V, M) is well defined. The kernel of this operator, called the heat kernel, 
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is defined by the equation 

(d t + F)U{t\x,x') = (21) 

with the initial condition 

U(0\x,xf) = 8(x,xf), (22) 

where 5(x, x') is the covariant Dirac distribution, the boundary condition 

Bi/;[U(t\x,x')]=Q. (23) 

and the self-adjointness condition 

U(t\x,x) =U*(t\x',x). (24) 

Hereafter all differential operators as well as the boundary data map act on 
the first argument of the heat kernel, unless otherwise stated. 

Let A be a complex number with a sufficiently large negative real part, 
Re A << 0. The resolvent can then be defined by the Laplace transform 

oo 

G(X) = [ dte tx U(t), (25) 



and by analytical continuation elsewhere. The heat kernel can be expressed, 
in turn, in terms of the resolvent by the inverse Laplace transform 



w+ioc 

1 



U(t) = — / dX e G(X) , (26) 



W — lOO 



where w is a sufficiently large negative real number, w << 0. As it has been 
done here, we will sometimes omit the space arguments if it does not cause 
any confusion. 



It is well known |3(| that the heat kernel U(t\x, x') is a smooth function 
near diagonal of M x M, i.e. for x close to x' , and has a well defined diagonal 
value 

U diag (t\x) = U{t\x,x), (27) 

and the functional trace 

Tr L2 exp(-tF) =J tr y f/ diag (t) , (28) 

M 
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where try is the fiber trace and the integration is defined with the help of 
the usual Riemannian volume element dvo\ g . 

It is also well known that in the smooth category the trace of the heat 



kernel has an asymptotic expansion as t — > + of the form [30] 



Tr L 2 exp(-tF) ~ t (k ~ m)/2 A k . (29) 

k=0 

Here are the famous so-called (global) heat-kernel coefficients (sometimes 
called also Minakshisundaram-Plejel coefficients). They have the following 
general form [ 3"0" |: 

M k = J 4? + J 41 (30) 

M dM 

A 2 k+i = J ( 31 ) 

where ajj, and a^, are the (local) interior and boundary heat-kernel coeffi- 
cients. The local interior coefficients a£ are also called HMDS (Hadamard- 
Minackshisundaram-De Witt-Seeley) coefficients in the literature. Hereafter 
the integration over the boundary is defined with the help of the usual Rie- 
mannian volume element d vol g on dM with the help of the induced metric 
9- 

The interior coefficients ajj. do not depend on the boundary conditions B. 
The even order coefficients calculated for Laplace-type operators up 

to [|I], [43| . The boundary coefficients do depend on both the operator 
F and the boundary operator B. They are far more complicated because 
in addition to the geometry of the manifold M they depend essentially on 
the geometry of the boundary dM. For Laplace-type operators they are 
known for the usual boundary conditions (Dirichlet, Neumann, or mixed 
version of them) up to Og [O, [T3], |3~3|. For oblique boundary conditions 



including tangential derivatives some coefficients were recently computed in 

However, the boundary value problem considered in the present paper 
with Zaremba type boundary operator B ( fL3|) is essentially singular. Even 
if the manifold M, its boundary dM and the operator F are all in smooth 
category, the coefficients of the boundary operator B are discontinuous on 
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Eo, which makes it a singular problem. For such problems the asymptotic 
expansion of the trace of the heat kernel has additional non-trivial logarithmic 



terms, \ T7[ |i~5|j 

oo oo 

Tr L 2exp(-tF B ) ~ ^ t^' 2 B k + log* ^ t fe/2 H k . (32) 

A:=0 k=0 

Whereas there are some results concerning the coefficients B k , almost nothing 
is known about the coefficients Hk. Since the Zaremba problem is local, or 
better to say 'pseudo-local', all these coefficients have the form 



2k 



4? + 



M 



B 



2k+l 



Si 



(,(!),! 
J 2k+1 



Hi, 



«,(!).! 
J 2k 



+ 



S 2 



(,(1),2 



l(2) 
^2fc ) 



S„ 



>(1),2 
y 2fc+l 



+ 



6 (2) 
u 2fe+l' 



So 



(33) 



(34) 



(35) 



Here the new feature is the appearance of the integrals over E , which com- 
plicates the problem even more, since the coefficients now depend on the 
geometry of the imbedding of the codimension 2 submanifold Eo in M that 
could be pretty complicated, even if smooth. The asymptotic expansion of 



the trace of the heat kernel has been studied recently in [ 40| . It has been 
shown there that the logarithmic terms do not appear, i.e. Hk = for 
any k, in the Zaremba type problem considered in the present paper, which 
confirmed the conjecture of ||. 



4 Parametrix of the Heat Equation 

Let us stress here that we are not going to provide a rigorous construction of 
the parametrix with all the estimates, which, for a singular boundary-value 
problem, is a task that would require a separate paper. For such a treatment 
the reader is referred to the papers |S8|, ESI E8J, [TIL |52| for the smooth 
case and to ]26|, [L7|, 0, |16|, [14], [15], [41], |37|, [27] for the singular case. 



Here we shall adhere instead to a pragmatic approach and will describe 
the construction of the parametrix that can be used to calculate explicitly 
the heat kernel coefficients B k as well as Hk- 
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4.1 Geometrical Framework 

First of all, we need to describe properly the geometry of the problem. Let 
us fix two small positive numbers £i,£2 > 0. We split the whole manifold in 
a disjoint union of four different parts: 

M = M int U M bnd = M int U M^ nd U M 2 bnd U M bnd . (36) 

Here Mq 1 " 1 is defined as the set of points in the narrow strip M bnd of the 
manifold M near the boundary dM of the width E\ that are at the same 
time in a narrow strip of the width e 2 near S 

M bnd = {x eM \ dist{x, dM) < e lt dist(x, S ) < e 2 } . (37) 

Further, M bnd is the part of the thin strip M bnd of the manifold M (of the 
width £i) near the boundary dM that is near Si but on the finite distance 
from So, i.e. 

M bnd = {x G M | dist(x, S x ) < £i, dist(x, S ) > e 2 } . (38) 

Similarly, 

M bnd = {x G M | dist(x, S 2 ) < ei, dist(x, S ) > £ 2 } • (39) 

Finally, M mt is the interior of the manifold M without a thin strip at the 
boundary dM, i.e. 

M int — M \ (M bnd U M bnd U M bnd ) = {x G M | dist(x, 9M) > £1} . (40) 

We will construct the parametrix on M by using different approximations 
in different domains. Strictly speaking, to glue them together in a smooth 
way one should use 'smooth characteristic functions' of different domains 
(partition of unity) and carry out all necessary estimates. What one has 
to control is the order of the remainder terms in the limit t — > and their 
dependence on e\ and e 2 - Since our task here is not to prove the form of the 
asymptotic expansion (|3"2"D, which is known, but rather to compute explicitly 
the coefficients of the asymptotic expansion, we will not worry about such 
subtle details. We will compute the asymptotic expansion as t — > in each 
domain and then take the limit £i,£2 — > 0. For a rigorous treatment see 
14] , |i~5| , 27| and the references therein. 
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We will use different local coordinates in different domains. In M mt we 
do not fix the local coordinates; our treatment will be manifestly covariant. 

In M\ nA we choose the local coordinates as follows. Let {e^}, (i = 
2, ...m), be the local frame for the tangent bundle TEi and x = (x l ) = 
(x 2 , . . . , x m ), (i = 2, ...,m), be the local coordinates on Ei. Let r = 
dist(x, Ei) be the normal distance to E 1; (r = being the defining equa- 
tion of Ei), and N = d r \ E be the inward pointing unit normal to E x . Then 
by using the geodesic flow we get the local frame {N(r, x), ej(r, x)} for the 
tangent bundle TM and the local coordinates x = (r,x) on M^ nd . The ge- 
ometry of Ei is described by the extrinsic curvature K (second fundamental 
form) 

W t e 3 = KijN , V t N = -K\ ej . (41) 

The coordinate r ranges from to £i, < r < £i. The local coordinates 
in M\ nd are chosen similarly. 

Finally, in M^ A we choose the local coordinates as follows. Let {e a (x)}, 
(a = 3, . . . , m), be a local frame for the tangent bundle TE and let x = 
(x a ) = (x 3 , . . . , x rn ) be the local coordinates on E . To avoid misunderstand- 
ing we should stress here that now we use the same notation x to denote 
coordinates on E (not on the whole DM). Let dist qm (x, E ) be the dis- 
tance from a point x on DM to So along the boundary dM. Then define 
y = + dist 8m(x, S ) > if x G Ei and y = — dist 3m(x, S ) < if x G E 2 . 
In other words, y = on E , (r = y = being the defining equations of 
E ), y > on Ei and y < on E 2 . Let n(x) = d y \ E be the unit nor- 
mal to E pointing inside Ei. Then by using the tangential geodesic flow 
along the boundary (that is normal to So) we first get the local orthonor- 
mal frame {n(y, x), e a (y, x)} for the tangent bundle TdM. Further, let the 
unit normal vector field to the boundary N(y, x) be defined as above. Then 
by using the normal geodesic flow to the boundary we get the local frame 
{N(r, y, x),n(r, y, x), e a (r, y, x)} for the tangent bundle TM and local coor- 
dinates (r,y,x) on MQ Dd . The geometry of E (codimension 2 manifold) is 
described by two extrinsic curvatures K and L and an additional vector T: 

V a e b = K ab n + L ab N . (42) 

V a n = -K\e h + T a N , V a N = -L b a e b - T a n . (43) 

The ranges of the coordinates r and y are: < r < e\ and —e^ < y < 
Finally, we introduce the polar coordinates 

r = p cos 9, y = p sin 9 . (44) 
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The angle 9 ranges from —tt/2 to tt/2 with 9 = —tt/2 on Ei and 9 = tt/2 on 
Eo. To cover the whole M bnd , p should range from to some £3 (depending 
on E\ and £ 2 ), < p < £3. 



4.2 Interior Parametrix 

This is the easiest case. The construction of the parametrix goes along the 
same lines as for manifolds without boundary (see, e.g. |19], [K], [j], |9|, f|. The 
basic case (when the coefficients of the operator F are frozen at a point xq) 
is, in fact, zero- dimensional, i.e. algebraic. By using the normal coordinates 
at Xq and Fourier transform one easily obtains the leading heat kernel 

K\t\x,x') = (47rf)- m / 2 exp (-^^) ■ (45) 

We try to find the fundamental solution of the heat equation near diagonal 
for small t, i.e. x — > x' and t — > + , that, instead of the boundary conditions 
satisfies asymptotic condition at infinity. This means that effectively one 
introduces a small expansion parameter e reflecting the fact that the points 
x and x' are close to each other and the parameter t is small. This can be 
done by fixing a point x$ = x' in M int , choosing the normal coordinates at 
this point (with g^ u (x') = S^ u ), scaling 

x — > x + six — x'), y — > x' + e(y — x'), t — > e 2 t, (46) 

and expanding in a power series in e. We will label the scaled objects by e, 
e.g. U £ . The scaling parameter e will be considered as a small parameter 
in the theory and we will use it to expand everything in power (asymptotic) 
series in e. At the very end of calculations we set e — 1. The non-scaled 
objects, i.e. those with e — 1, will not have the label e. Another way of 
doing this is by saying that we will expand all quantities in the homogeneous 
functions of (x — x'), (y — y') and \ft. This construction is standard and we 
do not repeat it here. 

One can also use instead a manifestly covariant method Jl9|, ||, [j], [3], f|, |43| , 
which gives a convenient formula for the asymptotics as t — > + 

00 

U int (t) ~ exp (-£) A 1/2 t (k - m)/2 a k , (47) 

k=0 
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where a = a(x,x') = (1/2) [ dist (x, x')] 2 is one half of the square of the 
geodesic distance between x and x', A = A(x,x') = g^ 1 ^ 2 {x)g^ 1 ^ 2 {x') 
det [—d*d* a(x, x')] is the corresponding Van Vleck-Morette determinant, 
g = detg^, and a/- = a>k(x,x') are the off-diagonal heat-kernel coefficients 
(note that odd order coefficients vanish identically, i.e. Ch2k+i = 0). These 
coefficients satisfy certain differential recursion relations which can be solved 
in form of a covariant Taylor series near diagonal [|TJ. 

The asymptotic expansion of the heat kernel on the diagonal reads 

oo 

^ g w-E t( ^ m)/2a f as ' ( 48 ) 

k=0 

where a^ iag (x) = ak{x,x). This asymptotic expansion can be integrated over 
the interior of the manifold M int . Since both the local interior coefficients 
and the volume element dvo\ g are regular at the boundary, these integrals 
have well defined limits as E\ — > 

diag / , diag 
'/, J 

M int M 



lim / tr v aT g = / tr v a™* . (49) 

£1^0+ 



Thus we obtain the local interior contribution to the global heat kernel co- 
efficients Bk. 

^ = trv4 aS - (50) 

As we already noted above all odd order coefficients vanish, fe^fe+i = 0- 
The explicit formulas for even order coefficients are known up to 
H], [43| . The first two coefficients have the well known form 

= (Air)- m/2 dimV , (51) 

bf = {A^)- m l 2 tv v [Q- l -R^ , (52) 
where R is the scalar curvature. 



4.3 Dirichlet Parametrix 

In this section we will follow closely the ideas of the paper . For an elliptic 
boundary- value problem the diagonal of the parametrix U h ^(t) in Mj 3Ild has 
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exponentially small terms, i.e. of order ~ exp(— r 2 /t), (recall that r is the 
normal geodesic distance to the boundary) as t — > + and r > 0. These terms 
do not contribute to the asymptotic expansion of the heat-kernel diagonal 
outside the boundary as t — > + . However, they behave like distributions 
near the boundary, and, therefore, the integrals over M^ nd , more precisely, 

£1 

the integrals lim £1 _ >0 / / dr(. . .), do contribute to the asymptotic expansion 

Si 

with coefficients being the integrals over Si. It is this phenomenon that 
leads to the boundary terms in the heat kernel coefficients. Thus, such terms 
determine the local boundary contributions to the global heat-kernel 
coefficients B^. The same applies to the Neumann parametrix and S 2 . 

The Dirichlet parametrix U hnd '^(t\x,x') in M] 3nd is constructed as fol- 
lows. Now we want to find the fundamental solution of the heat equation 
near diagonal, i.e. for x — > x' and for small t — > in the region M- bnd close 
to the boundary, i.e. for small r and r', that satisfies Dirichlet boundary 
conditions on Si and asymptotic condition at infinity. We fix a point on the 
boundary, x G S 1; and choose normal coordinates on S x at this point (with 
9ij(0,x ) = Sij). 

The basic case here (when the coefficients of the operator F are frozen 
at the point xq is one- dimensional. The zeroth-order term UQ nd '^ is defined 



by the heat equation 
where 



(d t + F )U^ (1) =0, (53) 
F = -d 2 r - d 2 , (54) 



the initial condition 

Uu nd ' {1) (0\r,x;r',x') = 5(r - r')5(x, x'), (55) 
the boundary conditions, 



u 



bnd,(l) 



= 0, (56) 

Si 



and the asymptotic condition 

lim f/ bnd ' (1) (t|r, x; r' , x') = lim f/ bnd ' (1) (t|r, x; r', x') = . (57) 



Note that the restriction to the boundary (...) applies only to the first 
argument, i.e. r — > 0. The operator F is a partial differential operator 
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with constant coefficients. By using the Fourier transform in the boundary 
coordinates (x — x') it reduces to an ordinary differential operator of second 
order. Clearly, the Eo part factorizes and the solution to the remaining one- 
dimensional problem can be easily obtained by using the Laplace transform, 
for example. The leading order Dirichlet parametrix thus has the form 

f/ bnd ' (1) (t\r, x; r', x) = K(t\r, x; r', x) - K(t\r, x; -r', x) (58) 

where 

- ' 1 u V -^J , (59) 

Note that in addition to the usual symmetry of the heat kernel, the Dirichlet 
parametrix possesses the following 'mirror symmetry 1 

U^ (1 \t\r,x;r',x') 

= -U^\t\ -r,x;r',x>) = -U^\t\r, x; -r', x') , (60) 

i.e. it is an odd function of the coordinates r and r' separately. 

To construct the whole parametrix, we again scale the coordinates. But 
now we include the coordinates r and r' in the scaling 

x — > x + e(x — xo), x — ^ xo + e(x — x ) (61) 

r — > er, r' — > er' , t — > e 2 t . (62) 
The corresponding differential operators are scaled by 

d -> -d, d r -> -d r , d t -»• \dt ■ (63) 

see 1 

Then, we expand the scaled operator F £ in the power series in e, i.e. 

oo 

F^F e ~Y,e»-*F n , (64) 

n=0 

where F n are second-order differential operators with homogeneous symbols. 
Since the Dirichlet boundary operator does not contain any derivatives and 
has constant coefficients on Si it does not scale at all. 
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The subsequent strategy is rather simple. We expand the scaled heat 
kernel in e 

oo 

[/bnd,(l) ^ ^ £ 2- m +r tf/ bnd,(l) ! ^ 
n=0 

and substitute into the scaled version of the heat equation and the Dirichlet 
boundary condition on Si. Then, by equating the like powers in e one gets 
an infinite set of recursive differential equations 

(d t + F )f4 bnd ' (1) = -J2 F n Utt {1 \ k = 1, 2, . . . , (66) 

n=l 

with the boundary conditions 

U^ {1 \t\0,x;r',x') = U^ (1 \t\r, x; 0, x') = 0, (67) 
and the asymptotic conditions 

lim U^ d ' (1 \t\r, x- r', x 1 ) = lim U^ d ' {1 \t\r, x; r', x') = . (68) 

In other words, we decompose the parametrix into the homogeneous parts 
with respect to (x — x ), (x' — x ), r, r' and \ft, i.e. 

U^ d ' {1 \t\r,x;r',x') 

= t {k-m)/2 U bod,(i) x' + r 1 / 2 (x-x'y,t- 1 / 2 r',x') , (69) 

in particular, on the diagonal we have 

U^ d ' (1 \t\r,x;r,x) = ^-^2^,(1) ^-^^-1/2^^ ? (7Q) 

and, therefore, 

00 

^dtg ,(1) w ~ X]^" m)/2c/ " nd ' (1) ( 1 l rl/2r ^; rl/2r ^) • ( 71 ) 

fc=0 

To compute the contribution to the asymptotic expansion of the trace of the 
heat kernel, we will need to compute the integral of ^diag'^W over M± nd . 
One should stress that the volume element should also be scaled 

OO L, 

dvol(r,x) — > dvo\(er,x) = dvo\(0,x) ■ /J ^77 (72) 

rv ! 

fc=0 
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where 

9k{x) = [d^F) \ ^ (73) 
Combining the above equations and changing the variable r = we obtain 

bnd,(i) f ,_ f ? , dvo\(r,x) bndi(1) 







E t(k ' m)/2 IJ2^9n(x) ! d^cu^-lm^,^), (74) 

fc=0 £ n=0 • ^ 

We note that even if the coefficients [/^ nd '^ 1 - ) satisfy the asymptotic regu- 
larity condition at r — ► oo (^8|) off-diagonal, the diagonal values of them do 
not fall off at infinity They have the following general form 

t4 bnd,(1) (iM;£,x) = m,*) + n (1) (£^) , (75) 

where Pk(£,,x) are polynomials in £ and : (£, x) are exponentially small, 
more precisely ~ £°exp(— £ 2 ) with some a, as £ — >• oo (which corresponds to 
t -> 0). 

Obviously, the integrals over the polynomial part over M 1 bnd vanish after 
taking the asymptotic expansion as t — > and the limit £i,£2 — > 0. The 
coefficients constitute simply the 'interior part' of the parametrix and are 
not essential in computing the boundary contribution. The coefficients K fc , 
in contrary, behave like distributions near Si. They give the Si contributions 
to the boundary heat kernel coefficients In the limit t — > the integral 

r 00 
Jo • •) becomes J cf£ (. . .) plus an exponentially small remainder term. 

o 

Then in the limit e\ — > we obtain integrals over Si up to an exponentially 
small function that we are not interested in. 

As the result we get the coefficients b% 1 in the form 

fc-i 00 

b k )A = E V / ^rtrv^i^x) . (76) 

n=0 n - I 
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These are the standard boundary heat kernel coefficients for Dirichlet bound- 
ary conditions. They are listed for example in |12|, [13] up to k = 4. The first 
two have the form 



0, 



6p = _( 4vr )-(™-i)/2 dimV^ 



( 4vr )-"V2 fomV-K 
3 



(77) 



where K is the trace of the extrinsic curvature (second fundamental form) 
of the boundary. 



4.4 Neumann Parametrix 

The construction of the Neumann parametrix in M\ nd is essentially the same 
except that now the boundary operator, in fact the endomorphism A, is not 
constant and should be also scaled, so that the scaled boundary conditions 
are §, g 

'1 



-d r + A £ ) Lp 



0, (78) 



where 

oo 

A^J^A*. (79) 

fe=0 

The zeroth-order operator Fq is given by the same formula ([54]) and the 
zero order boundary operator is just the standard Neumann one. The basic 
zero-order problem can again be easily solved by 

£ybnd,(2) ^ ^ _ ^£| r? £. r >^ _|_ K(t\r, x; — r', x) (80) 



with the same kernel K d59|). Note that the Neumann parametrix has another 
mirror symmetry 

U^ d ' (2) (t\r,x;r',x') 

= (t\ _ r> x - r ' ? £') = u^P) ( t \r, x; -r', x') , (81) 

i.e. it is an even function of the coordinates r and r' separately. 
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The construction of the parametrix goes along the same lines as in Dirich- 
let case. We have the recursive differential equations 

(d t + F )U^ d > {2) = -J2 FnUtT\ k — 1,2,..., (82) 



n=l 



with the boundary conditions 



fc-i 



d r U k 



bnd,(2) 



n=l 



-E A ^ b - d i 2 i v ' ( 83 ) 



and the asymptotic conditions 

lim £/knd,(2) ^. ^ _~/^ _ jjbnd,(2) ^. ^ ^/^ _ g ^ ( S4 ) 

r — >oo r 1 — >oo 

As we already noted above the restriction to the boundary applies only to 
the first argument r. One can repeat here everything said at the end of the 
previous subsection about Dirichlet parametrix. We have again homogeneity 
property 

U^ 2 \t\r,x;r',x') 

= )/ 2 c/ fe bnd ' (2) (i|r 1 / 2 r, x' + r x i\x - x'); t^'V , x') (85) 

and the following expansion for the diagonal 

oo 

^ag (2 V) ~ ^* (fc " m)/2c/ if nd,(2) ( 1 l rl/2r ' £ ; rl/2r ' £ ) • ( 86 ) 

A:=0 

By separating the polynomial and exponentially small parts, 

U^ 2 \l\^x^x) = Pk{^x) + Yi 2 \£,x) , (87) 

and repeating the arguments at the end of the previous subsection we obtain 
the S 2 contributions to the boundary heat kernel coefficients b k 

k-l 00 

^ = E ~y9n [ «r try^iUtf,*) . (88) 
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These are the standard boundary heat kernel coefficients for Neumann bound- 
ary conditions. They are listed for example in |12|, [13] up to k = 4. The first 
two have the form 



■,(1).2 







b^ 2 = (4vr)^ m - 1 )/ 2 dim\/^ 



l(1)>2 



Un)- m/2 dim V-K. 
v ; 3 



(89) 



4.5 Mixed Parametrix 

This is the most complicated (and the most interesting) case, since here 
the basic problem with frozen coefficients on E is two-dimensional. More 
precisely, in MQ nd the basic problem is on the half-plane. Since the origin is 
a singular point, we will work in polar coordinates introduced above. 



4.5.1 Basic Problem (Zeroth Order) 



First of all, we need to solve the basic problem for operators with frozen 
coefficients at a point xq on Eo. We choose normal coordinates on Eo at this 
point (with g a b(0, 9, xq) = S a b) and the polar coordinates (p, 6) in the normal 
bundle described above. Then the zero order operator F$ has the form 



~d P - - 2 d 2 - d 2 
p p 2 



(90) 



where d 2 = g ab d a db. The zero order inward pointing normal N to the bound- 
ary in polar coordinates has the form 



N 



0, 



y>0 p 



p>0, 



(91) 



N 



T.2 



0, 



y<0 p 



p>0, 



(92) 



iVn 



So 



Or 



y=0 



d,. 



p=o, e=o 



(93) 



Now the boundary operator is discontinuous, and there is a singularity at 
the origin p = 0. 
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Again the part due to S factorizes 
U^°\t\p,9,x;p',9',x') 
exp 



(47rt) -(m-2)/2, 



X — X 



r\2 



At 



V(t\p,9;p',9'), 



(94) 



where ^(t\p,9; p' ,9') is a two-dimensional heat kernel. It is determined by 
the heat equation 



(dt - d 2 p - -d p - y 2 d 2 ^ 9(t\p, 9; p', 9') = , 



the initial condition 

tt(O + |p,0;p',0' 
the boundary conditions 



'pp 



_ 5 (p - p ')5(9 - 9') , 



V(t\p,9;p',9') 
d e *(t\p,9;p',9') 



0, 



= 0, 



the symmetry condition 

y(t\p,9;p',9') = y(t\p',9';p,9), 

as well as some boundary conditions at p — > + and p — > oo. 
We require certain regularity conditions at infinity 



oo 

J dp^p~p~' \^(t\p,9;p',9')\ <oo, 



in particular, 



lim ^/pp'^(t\p,9;p',9') = lim d p y/pp* V(t\p,9; p' ,9') 



0. 



(95) 

(96) 

(97) 
(98) 

(99) 



(100) 



(101) 



As far as the boundary condition on S , i.e. at p = 0, is concerned, 
we will see that it cannot be a generic condition, rather, like for the usual 
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Frobenius theory of differential equations near singular points, one has a 
couple of possibilities for the type of singularity that need to be specified. 

Since we are looking for the solution of the heat equation whose diagonal 
is integrable near boundary, we require 



dpp\^(t\p, 9;p, 0)| < oo 



(102) 



This effectively imposes a restriction on the type of the singularity at p — > 0, 
i.e. the singularity of the heat kernel at p — > must be weaker than (pp')^ 1 . 
More precisely, we assume that 



^\y(t\p,6;p',6')\<oo. 



(103) 



We will see that this is still not enough to fix a unique solution and 
another boundary condition at p — > is needed. Since the point p = is 
singular, this boundary condition cannot be imposed arbitrarily. Also it does 
not follow from the boundary conditions on Ei and E 2 . We impose it in one 
of the following forms 



^*(t|p,0;p',0') 







or 



o 



p=ch 



(104) 



(105) 



where s is a real parameter. We will see that the heat kernel asymptotics do 
depend on this boundary condition as well. We call the boundary condition 
(|104j) "regular" boundary condition. It corresponds formally to the limit 
s — > oo. 

To construct the heat kernel we study first the operator 

L = -dl 



on the interval [— f , f ] with the boundary conditions 



<p{0) 



0. 



0. 



(106) 



(107) 



It is not difficult to find the spectral resolution of this operator. Its orthonor- 
mal eigenfunctions and eigenvalues are 



<Pn(9) 



COS 



n + 



1 







71 



108) 



An = [n + 



where n — 0, 1, 2, 

By separating the variables 

oo 

*(t\p,9;p',9') = Y,Vn{6) Vn {6')u n (t\ P -,p' 
we obtain the equation 



71=0 



P P 2 



with the initial condition 



1 



Ur X0 + \p;p') = —5(p-p'), 
the symmetry condition 

u n (t\p;p') = u n (t\p';p) 
and the asymptotic conditions at infinity 



oo 

J dp^fpp 1 \u n (t\p,9;p',9')\ < oo, 



lim yfpf) u n (t\p; p) = lim d p yfpp 1 u n (t\p; p') 

p— >oo p^oo L 

The boundary conditions at p = are 
and, more precisely, one of the following 



or 



^u n (t\p,9;p',9') 

- p=0+ 

(dp-s) y/pp 1 ; u n (t\p,9;p',9') 



= 0, 



p=0+ 
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Let us consider the operator 



1„ If 1 N 2 



D n = -d 2 p --d P + y 2 i^ + -j . (119) 

It has the "eigenf unctions" J u {pp): 

D n J u {pp) = p 2 J u (pp) , (120) 

where p is a positive real parameter, J„(z) are Bessel functions of the first 
kind of order u, and v can take one of two values, either v — (n + 1) or 
v = — + However, the behavior at p — ► of the Bessel function 

J( n+ i\{pp) for ri > 1 is too singular, ~ p~( n+1 z\ which violates the inte- 

grability condition ( |116| ). This means that for any n > 1 we have to choose 



v — (n + I) . Note that these are not "true" eigenfunctions, since they are 
non-normalizable. Rather they satisfy the following "orthogonality" condi- 
tion 

oo 

/ dppJ n+ i(pp)J n+ i(pp) = —=5(p- p) . (121) 
o 

In the case n = both choices are possible, i.e. v = +1/2 or v = —1/2, 
which makes the analysis of the problem more complicated. Therefore, we 
will treat the cases n > 1 and n = separately. 

Case I. We consider first the case n > 1. We will solve this problem by 
employing the Hankel transform which is well defined in the class of functions 
satisfying the conditions imposed above. We define 

oo 

v n (t\p,p) = J dp pj n+ i(pp)u n (t\p, p') , (122) 
o 

then 

oo 

u n (t\p,p')= / dpp J n+ i(pp)v n (t\p, p') , (123) 







Next, by integrating by parts and using the eq. ( |120| ), we compute the 
Hankel transform 

oo oo 
J d PpJn+\{Pp) D nU n {t\p,p')=p 2 J dp p J n+ i(fip)u n (t\p, p) 
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+ P 



d p J n Mpp) u n (t\p,p') - J n+ i(fip)d p u n (t\p, p') 



(124) 



Finally, by taking into account the boundary conditions ( |115[ ) and ( |116[ ), and 
the asymptotic form of the Bessel functions, we obtain 



(125) 



dppJ n+ i(pp)D n u n (t\p,p') = p 2 v n (t\p,p') 
Thus, the Hankel transform of the heat equation (|111 ) is 

{d t + p 2 )v n {t\p, P ') = o. 



From (|112|) we also obtain the initial condition 

v n (0 + \p,p') = J n+ i{pp) . 
It immediately follows that 

v n (t\p,p') = e'^ 2 J n+ i(pp') 

and, therefore, 



(126) 



(127) 



(128) 



oo 

u n (t\p,p') = J dp jue"*" 2 J n+1 / 2 (pp)J n+1/2 {pp') ■ (129) 



This integral can be computed by using the properties of the Bessel functions. 
We obtain finally 



u n {t\p, p) = ^exp 



p 2 + p 12 

At 



In+l/2 



pp. 

2t 



(130) 



where I n+ i/2(z) is the modified Bessel function of first kind. Note that this 
solution satisfies both boundary conditions (|117|) and ( |118|) . 

Case II. Now let us consider the case n = 0. As we have seen the condi- 
tion of integrability (|116|) does not fix the solution uniquely, since there are 
two linearly independent solutions that satisfy that condition, which corre- 
sponds to the choices v — — 1/2 and v = +1/2. The Hankel transform in this 
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case reduces to the standard cosine and sine Fourier transforms. However, 
we will not use them, but will solve the heat equation directly. 
Let us single out the allowed singular factor 



u (t\p,p') 



=w(t\p, p') . 



(131) 



Then, the heat equation ( |1 1 1| ) , the initial condition ( |112| ), and the boundary 
conditions (|117|) and ( |118| ) take the form 



(d t - d 2 p )w(t\p, p') = , 
w(0 + |p,p) = 5{p- p) 



w(t\p,p) 



p=0 



or 



(dp - s)w(t\p,p) 



0. 



p=0 







(132) 
(133) 
(134) 

(135) 



There is also the usual regularity condition at infinity p — > oo. 

As we see, w is just the standard one-dimensional heat kernel on the half- 
axis. By using the Laplace transform we easily obtain the solution of this 
problem 



w(t\p,p') 



1 

2vri 



c+ioo 



d\ e 



-tx 



2 V /Z A 



exp 



-\f-X\p- p'\ 



-A 



-A + s 



exp 



-v^A(p + p') 



(136) 



where c is a sufficiently large negative real constant, i.e. y— c > — s, and 
y/— A is defined in the complex plane of A with a cut along the real positive 
half-axis, so that Re y/—\ > 0. Notice that the boundary conditions (|134|) 
correspond to the limit s — > +oo. The limit s — > — oo is not well defined 
since the constant c depends on s and would have to go to — oo as well. 
Next, let us change the variable A according to 



A 



(137) 



where Imp > 0. In the upper half-plane, Imp > 0, this change of variables 
is single-valued and well defined. Under this change the complex A-plane is 
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mapped onto the upper half /i-plane, and the cut in the complex A-plane 
along the positive real axis from to oo is mapped onto the whole real axis 
in the /i-plane. 

The contour of integration in the complex //-plane is a hyperbola going 
from (e* 37r,/4 )oo through the point y/—c to (e l7r / 4 )oo. It can be deformed to a 
contour C that is above all poles of the integrand. It comes from — oo along 
the real axis, encircles posible poles on the imaginary axis in the clockwise 
direction, and goes to +00 along the real axis. 

After such a transformation we obtain 



w(t\p, p') 



S 2n { 6XP + ~ 

c 

exp [— tp 2 + ip(p + p' 



, , (138) 

p + ts L v ' J J v ' 

This function is an analytic function of s since the contour C is above the 
pole at —is. Therefore, we can compute it, for example, for s > 0, and then 
make an analytical continuation on the whole complex s-plane. So, let s > 0. 
Then the pole —is is in the lower half-plane. Therefore, the contour C can 
be deformed to just the real axis, i.e. —00 < p < 00. Next, we use the 
following trick 



^ = 1-2*^- 

p + is p + is 



1 — 2s / dpe 



(139) 



This integral converges since s > 0. Substituting this equation in ( |138| ) and 
evaluating the Gaussian integral over p, we obtain 

{p-m , r (p+pT 



w{t\p,p') 



(Ant) 



-1/2 



exp 



2s / dp exp 



At 

jp + p'+p) 2 
At 



+ exp 



At 



ps 



(140) 



which can be expressed in terms of the complimentary error function 

(p-p') 2 i . r (p+pT 



w(t\p, p') = (Ant)- 1 ? 2 



exp 



At 



+ exp 



At 



-2y/7TsViexp [ts 2 + (p + p')s] erfc + 
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Here erfc (z) is defined by 



oo 

2 f 2 
erfc (z) = —j= / du e~ u . 

V n J 



(142) 



The case s < can be analyzied either directly or by the analytical contin- 
uation. The direct computation is different since now the pole —is is in the 
upper half-plane and one has to take into account the residue at this pole. 
However, the integral along the real axis is also different, so that the sum is 
the same. In other words, the result for s < has the same analytical form 

Finally, we obtain the heat kernel component Uq: 



u (t\p,p') = (Ant)- 1 / 2 - 



'PP 



exp 



(p ~ P'f 



to 



+ exp 



(p + P 



A2 



to 



— 2\/vrsV^exp [ts 2 + (p + p')s] erfc 
In the particular case s = we get 



P + P' 
2Vi 



+ sVi 



(143) 



u (t\p,p') = (47rf) 



-1/2. 



'PP 



exp 



t\2 



(P ~ P ') 

4t 



+ exp 



(P + P 



l\2 



to 



(144) 

The case s — > +oo corresponds to the regular boundary conditions ( 117|) . In 
this case the solution reads 



U (t\p, P ') = (47Tt)- 1 / 2 - 



1 



PP 



exp 



{P-P 



l\2 



2t 



exp 



p 2 + p 12 

to 



1 1/2 



it 
P± 

2t 



exp 



(P + P 



A2 



to 



(145) 

obtained by the Hankel 



and coincides with the solution ( |130| ) for n 
transform. 

Combining our results and using the explicit form of the eigenfunctions 
(p n , we obtain the heat kernel 



^(t\p,9;p',9') = (Ant)- 1 mp,p') 
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X < cos 



9-9' 



+ (Ant) exp 



x )n[^,9-9')+n[^,9 + 9' + 7i 



+ cos 



p 2 + p' 2 



9 + 9' + 7T 



At 



2t 



PP' 



2t 



(146) 



where 





r ip+p'f] 


(£)' { exp 


At 





-^/nsViexp [ts 2 + (p + //>] erfc (^/=- + S ^J J > ( 147 ) 

r 



and 



^0,7) = 2^/ n+ i /2 (z)cos + 

n=0 L V / 



7 



(148) 



Notice that for the "regular" boundary conditions ( |104| ), which correspond 
to the limit s — > +oo, the function p') vanishes. 

This series can be evaluated by using the following integral representation 
of the Bessel function 



4+1/2(2) 



dpe pz (l - p 



2\n 



V2 

Substituting this integral in the series and summing over n we obtain 



(149) 



n(*, 7 ) 




2n 



dpe pz < exp 



\{1-P 2 )ze^ + U 1 



-1 



+ exp 



i(l-p=»)ze-*r-Iry 



(150) 



The remaining integral can be expressed in terms of the error function, so 
that finally we get 



fi(3, 7 ) = e 2COS7 erf 



1 : C( »s ( — 



(151) 
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where the error function is defined by 

2 

erf (z) = 4= / dpe~ p2 . (152) 

V 71 J 
o 

By adding the S factor we obtain the final answer for the parametrix 
f/ bnd ' (0) (t|p, 0; p', 9') = L(t\p, 0, x; p', 9', x') + L(t\p, 0, x; p' , -9' - vr, x') (153) 
where 

L(t\p,9,x;p',9',x') 

= (47rt)-™/ 2 exp (-^^) *(%p')cos 

+ (4vrt)- m / 2 exp j-1 [|x - x'\ 2 + p 2 + p' 2 - 2pp' cos(0 - 9')] j 
X(T f I 7^ C0S (^zZ ) j . (151) 
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An important corollary from this formula are the symmetries of the heat 
kernel. First of all, we have the usual 'self-adjointness' symmetry 

0->0', P^P- (155) 

Second, we have the 'periodicity' symmetries 

0^0 + 4vrn, 0'->0' + 47rm, n,meZ. (156) 

Finally, there is additional 'mirror' symmetry 

0^0, 0' -> -0' - 7T (157) 

9 -> -0 - vr, 0' -> 0' (158) 

Note the essential difference of the symmetries of the mixed parametrix versus 
those of the Dirichlet and Neumann parametrices. The mixed parametrix is 
a periodic function of the angles (expected), but not with the period 27T but 
with the period 4-7T (not expected). That is why there are two different mirror 
images, (p, —9 — n,x) and (p, —9 + ir, x), of a point with the coordinates 
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(p,9,x). In other words the double reflection of a point does not bring it 
back — the double image is not identical with the original point. Denoting by 
T the transformation 9 — > —9 — it we have 



T 4 =Id, but TV W. (159) 

The operator T has four eigenvalues 1, —1, % and —%. Whereas the first two, 
1 and —1, are the standard ones, the latter two, % and —i, correspond to 
some new images. This might have some interesting applications. 
The diagonal of the mixed parametrix is easily found to be 



/-rbnd,(0) / ,x 
U diag,0 \ L ) 



exp 



+(l-sin0) 



(Anty m/2 

p 2 cos 2 9 
t 

4 Vt 



erf 



erfc ( -^= 



psin 9 



A p 



exp 



-y/nsVi exp (ts 2 + 2p s) erfc ^-^= + sVi^j | . 



(160) 



Now we compute the integral of the diagonal of the parametrix over Mq nd 

7I-/2 



£3 



try*/. 



bnd,(0) 



diag,0 



(*) 



M bnd 



dpp J Mtr v u£2™(t) 

-it/2 



(161) 



for some finite e 3 > \Je\ + e\ > 0. 

First of all, obviously the integrals over 9 of the odd functions in 9 vanish 
identically. So, we only need to consider the even part. Second, since in the 
limit £3—^0 the volume of Mq nd vanishes, the regular part of the heat kernel 
diagonal does not contribute to the trace either. It is only the singular part 
of the heat kernel diagonal, which behaves like a distribution near S , that 
contributes to the integral in the limit £3^0. 

The integral over p can be computed exactly. It reads 

_2 



/ 



tr v U 



bnd,(0) 



diag,0 



= J(4nty 



M bnd 



+ t 



-^ + 2tt0 
4 



Vii 



X(t) 



(162) 
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where 

0(z) = e z2 erfc(z), (163) 

and 



X(t) = 2^3 ^(-fj + ^-^j erfc ^ 

-27rtexp (ts 2 + 2se 3 ) erfc + Vis^j . (164) 

Notice that vre|/2 is nothing but the area of the semi-circle of radius £3, so 
that vol (S )7re2/2 = vol(M bnd ). In the limit e 3 -> this term does not 
contribute to the asymptotics. 

By using the asymptotic behavior of the error function as z — > 00 

1 _ 2 

erfc (z) ~ —^e- z ■ (165) 

we find that the function X(t) is exponentially small, i.e. it is suppressed by 
the factor ~ exp(— £§/£), as t — > 0, and, therefore, does not contribute to the 
asymptotic expansion of the heat kernel in powers of t fl32|) either. 

The behavior of the function Q(\/is) depends on the parameter s. For a 
finite s in the limit t — > we have 

6 (y/ts\ = l+0(t 1/2 ). (166) 

It immediately follows that for a finite s, i.e. for the boundary conditions 
( |105|) the singular heat kernel coefficient b% is equal to 

b {2) = ( 4vr )-(™-2)/2 dim V ]_ ( 167 ) 

16 

Notice that it does not depend on s explicitly. 

Finally we analyze the regular boundary conditions ( |104j ), which corre- 
sponds formally to the limit s — > +00. By using ( |165|) we see that for a finite 
t the function 0(\/ts) vanishes in the limit s — > 00: 

©(v^s) =0. (168) 
Therefore, in this case the coefficient is 



= _(4vr)- (m - 2)/2 dim 7 — . (169) 



2 16 



(2) 

We see that the heat kernel coefficients b k depend on the type of the 
additional boundary conditions, i.e. ( |104| ) vs. ( |105|) , at S . 
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5 Conclusions 

We have studied Zaremba boundary value problem for second-order par- 
tial differential operators of Laplace type, when the manifold as well as the 
boundary are smooth and the differential operator has smooth coefficients 
but the boundary operator is discontinuous on the boundary, it jumps from 
the Dirichlet to Neumann type boundary operator. Since this problem is 
not smooth there could be additional logarithmic terms in the asymptotic 
expansion of the trace of the heat kernel (0). However, Seleey ||0] has 



shown recently that such terms do not appear and there is classical asymp- 
totic expansion in half-integer powers of t only. This seems to contradict the 



conclusions of pl| , where it has been shown that such an expansion with 
locally computable coefficients does not exist. The term 'locally computable' 
is confusing though. As we have seen the calculation of the coefficients of the 
asymptotic expansion of the trace of the heat kernel involves the knowledge 
of some global information, i.e. the spectrum ( |109| ) of the operator L (|106|) 
with mixed boundary conditions. So, one could say that these coefficients are 
locally computable in the coordinates x and p but are global in the coordinate 
9. Thus the standard asymptotic expansion in powers of t (without logarith- 
mic terms) still exists with coefficients (|3"3"D, ( |5^ ) given by integrals over M, 
Si, £2 and So- The interior coefficients, bf\ the co-dimension one coeffi- 
cients, typ' 1 and , are 'locally computable', but the co-dimension two 
(2) 

coefficients, b k , are 'global' in 9 (or pseudo-local) and require new methods 
of calculation (e.g. like the approach of this paper). They are constructed 
from the local invariants on £ . It is the numerical coefficients that are 
global. 

Let us formulate briefly our main results. First of all, we provide the 
correct formulation of the Zaremba type boundary value problem. We find 
that the boundary conditions on the open sets Si and £2 are not enough to 
fix the problem, and an additional boundary condition along the singular set 
E is needed. This additional boundary condition can be considered formally 
as an 'extension' of Dirichlet conditions from Si to S , or an 'extension' of 
Neumann conditions from E2 to Eo- However, strictly speaking the boundary 
conditions on Eo does not follow from the boundary conditions on E \ Eo and 
can be chosen rather arbitrarily. One needs some additional 'physical' crite- 
ria to fix this boundary condition. Second, we describe the geometry of the 
problem, which involves now some nontrivial geometrical quantities (normal 
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bundle and extrinsic curvatures) that characterize properly the imbedding 
of a co-dimension two submanifold S in M. The higher order coefficients 

(2) 

b k are invariants constructed from those geometric quantities. Next, we 
describe the construction of the parametrix of the heat equation in the in- 
terior of the manifold M, in a thin shell close to Si and S 2 , and finally, in 
a thin strip close to E . We used the standard scaling device; the differ- 
ence is just in what coordinates are involved in scaling. Finally, we have 
explicitly found the off-diagonal parametrix in Mq 11 * 1 , the thin strip close to 
S , in the leading approximation, and used it to compute the first nontrivial 
'global' coefficient, b 2 , of the heat kernel asymptotic expansion. We con- 
sidered two types of the additional boundary condition along S , one being 
the 'extension' Dirichlet boundary conditions (that we called regular bound- 
ary condition), and another being the 'extension' of the Neumann (or rather 
Robin) boundary conditions. We have shown that the result, i.e. the coeffi- 

(2) 

cient b 2 , does depend on the type of the boundary condition, i.e. Dirichlet 
vs Neumann, but does not depend on the parameter s of the Robin boundary 
condition (it will however contribute to the higher order coefficients). 

Acknowledgements 

I would like to thank Jochen Briining, Stuart Dowker, Giampiero Esposito, 
Stephen Fulling, Peter Gilkey, Gerd Grubb, and Werner Miiller for stimu- 
lating and fruitful discussions. I am also very grateful to Robert Seeley for 
clarifying discussions of the boundary conditions and sharing the preliminary 
results. The support by the NSF Block Travel Grant DMS-9988119, by the 
MSRI, and by the Istituto Italiano per gli Studi Filosofici and the Azienda 
Autonoma Soggiorno e Turismo, Napoli, is gratefully acknowledged. 

References 

[1] Avramidi, I. G.: A covariant technique for the calculation of the one- 
loop effective action, Nucl. Phys. B355 (1991), 712-754. Erratum: 
Nucl. Phys. B509 (1998), 557-558. 

[2] Avramidi, I. G.: A method for calculating the heat kernel for mani- 
folds with boundary, Yadernaya Fizika 56 (1993), 245-252, [Russian]; 
Physics of Atomic Nucleus 56 (1993), 138-142 [English] . 



37 



[3] Avramidi, I. G.: Green functions of higher-order differential operators, 
J. Math. Phys. 39 (1998), 2889-2909. 

[4] Avramidi, I.G.: Covariant techniques for computation of the heat ker- 
nel, Rev. Math. Phys. 11 (1999), 947-980. 

[5] Avramidi, I. G.: Heat kernel asymptotics of non-smooth boundary 
value problem, New Mexico Tech, 1999; in: M. van den Berg and 
V. Liskevich (eds), Workshop on Spectral Geometry, Abstracts of Int. 
Conf., University of Bristol, Bristol, UK, July 10-15, 2000. 

[6] Avramidi, I. G. and Esposito, G.: Lack of strong ellipticity in Euclidean 
quantum gravity, Class. Quant. Grav. 15, (1998), 1141-1152. 

[7] Avramidi, I. G. and Esposito, G: Gauge theories on manifolds with 
boundary, Comm. Math. Phys. 200 (1999), 495-543. 

[8] Avramidi, I. G. and Esposito, G.: Heat kernel asymptotics of the 
Gilkey-Smith boundary value problem, in: V. Alexiades and G. Siopsis 
(eds), Trends in Mathematical Physics, AMS/IP Studies in Advanced 
Mathematics, vol. 13, American Mathematical Society and Interna- 
tional Press, 1999, pp. 15-34. 

[9] Avramidi, I. G. and Schimming, R: Algorithms for the calculation of 
the heat kernel coefficients, in: M. Bordag (ed), 'Quantum Field Theory 
under the Influence of External Conditions ', Teubner-Texte zur Physik, 
vol. 30, Teubner, Stuttgart, 1996, pp. 150-162. 

[10] Berline, N., Getzler, E. and Vergne, M.: Heat Kernels and Dirac Op- 
erators, Springer- Verlag, Berlin-New York, 1992. 

[11] Booss-Bavnbek, B. and Wojciechowski, K. P.: Elliptic Boundary Prob- 
lems for Dirac Operators, Birkhauser, Boston, 1993. 

[12] Branson, T. and Gilkey, P. B.: The asymptotics of the Laplacian on a 
manifold with boundary, Comm. Part. Diff. Eq. 15 (1990), 245-272. 

[13] Branson, T. P., Gilkey, P. B., Kirsten, K. and Vassilevich, D. V.: Heat 
kernel asymptotics with mixed boundary conditions, Nucl. Phys. B563 
(1999), 603-626. 



38 



[14] Briining, J. and Seeley, R.: Regular singular asymptotics, Adv. in 
Math. 58 (1985), 133-148. 

[15] Briining, J. and Seeley, R. T.: The expansion of the resolvent near a 
singular stratum of conical type, J. Fund. Anal. 95 (1991), 255-290. 

[16] Callias, C: The heat equation with singular coefficients I, Comm. 
Math. Phys. 88 (1983), 357-385. 

[17] Cheeger, J.: On the spectral geometry of spaces with cone-like singu- 
larities, Proc. Nat. Acad. Sci. USA 76 (1979), 2103-2106. 

[18] Cheeger, J.: Spectral geometry of singular Riemannian spaces, J. Diff. 
Geom. 18 (1983), 575-657. 

[19] De Witt, B.S.: The Spacetime Approach to Quantum Field Theory, in: 
B. S. De Witt and R. Stora, R. (eds.) Relativity, Groups and Topology 
II, North Holland, Amsterdam, 1984, pp. 383-738. 

[20] Dowker, J. S.: The N U D problem, University of Manchester, 2000, 
|hep-th/0007i27| . 

[21] Dowker, J. S., Gilkey, P. B. and Kirsten, K.: On properties of the 
asymptotic expansion of the heat trace for the N/D problem, Int. J. 
Math. 12 (2001), 505-517. 

[22] Dowker, J. S. and Kirsten, K.: Heat-kernel coefficients for oblique 
boundary conditions, Class. Quant. Grav. 14 (1997), L169-L175. 

[23] Dowker, J. S. and Kirsten, K.: The 03/2 heat-kernel coefficient for 
oblique boundary conditions, Class. Quant. Grav. 16 (1999), 1917- 
1936. 

[24] Elizalde, E. and Vassilevich, D. V.: Heat Kernel Coefficients for Chern- 
Simons Boundary Conditions in QED, Class. Quant. Grav. 16 (1999), 
813-823. 

[25] Fabrikant, V. I.: Mixed boundary value problems of potential theory and 
their applications in engineering, Kluwer, Dordrecht, Boston, 1991. 



39 



[26] Fedosov, B. V.: Asymptotic formulas for the eigenvalues of the Laplace 
operator in the case of a polyhedron, Soviet Math. Dokl. 5 (1964), 988- 
990. 

[27] Gil, J. B.: Full asymptotic expansion of the heat trace for 
non-self-adjoint elliptic cone operators, Temple University (2001), 
|math.AP/0004161| . 

[28] Gilkey, P. B. and Smith, L.: The eta invariant for a class of elliptic 
boundary value problems, Comm. Pure Appl. Math. 36 (1983), 85-132. 

[29] Gilkey, P. B. and Smith, L.: The twisted index theorem for manifolds 
with boundary, J. Diff. Geom. 18 (1983), 393-344. 

[30] Gilkey, P. B.: Invariance Theory, the Heat Equation and the Atiyah- 
Singer Index Theorem, Chemical Rubber Company, Boca Raton, 1995. 

[31] Grubb, G.: Properties of normal boundary value problems for elliptic 
even-order systems, Annali della Scuola Norm. Sup. di Pisa, vol. 1, 
ser. IV (1974) 1-61. 

[32] Grubb, G.: Functional Calculus of Pseudodifferential Boundary Prob- 
lems, Progress of Mathematics, 65, Birkhauser, Boston, 1996. 

[33] Karol', A. I.: Asymptotics of the parabolic Green function for an ellip- 
tic operator on a manifold with conical points, Math. Notes 63 (1998), 
no. 1-2, 25-32. 

[34] Kirsten, K.: The heat kernel coefficient on a manifold with boundary, 
Class. Quant. Grav. 15 (1998), L5-L12. 

[35] Lesch, M.: Operators of Fuchs type, conical singularities, and asymp- 
totic methods, Teubner-Texte zur Math., vol. 136, B.G. Teubner, 
Stuttgart, Leipzig, 1997. 

[36] McAvity, D. M.. and and Osborn, H.: Asymptotic expansion of the 
heat kernel for generalized boundary conditions, Class. Quant. Grav. 
8 (1991), 1445-1454. 

[37] Mooers, E. A.: Heat kernel asymptotics on manifolds with conic sin- 
gularities, J. Anal. Math. 78 (1999), 1-36. 



40 



[38] Seeley, R. T.: Topics in pseudo-differential operators, in: CIME Con- 
ference on Pseudo-Differential Operators 1968, Edizioni Cremonese, 
Roma, (1969), 169-305. 

[39] Seeley, R. T.: The resolvent of an elliptic boundary value problem, Am. 
J. Math. 91 (1969), 963-983. 

[40] Seeley, R. T.: Trace Expansions for the Zaremba Problem, University 
of Massachusetts at Boston, 2001; submitted to Comm. Part. Diff. Eq. . 

[41] Simanca, S. R.: Mixed elliptic boundary value problems, Comm. Part. 
Diff. Eq. 12 (1987), 123-200. 

[42] Sneddon, I. N.: Mixed boundary value problems in potential theory, 
Wiley, New York, 1966. 

[43] van de Ven, A. E. M.: Index free heat kernel coefficients, Class. Quant. 
Grav. 15 (1998), 2311-2344. 



